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We investigate the thermodynamics of model systems exhibiting two-scale fractal spectra. In 
particular, we present both analytical and numerical studies on the temperature dependence of the 
vibrational and electronic specific heats. For phonons, and for bosons in general, we show that the 
average specific heat can be associated to the average (power law) density of states. The corrections 
to this average behavior are log-periodic oscillations which can be traced back to the self similarity 
of the spectral staircase. In the electronic case, even if the thermodynamical quantities exhibit 
a strong dependence on the particle number, regularities arise when special particle numbers are 
considered. Applications to substitutional and hierarchical structures are discussed. 
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I. INTRODUCTION 

In the last 15 years quasiperiodic structures have been 
studied intensively. Aside from their purely theoretical 
interest, these studies were in part motivated by the dis- 
covery of the quasicrystalline state JIJ together with the 
possibility of experimental realization of quasiperiodic 
superlattices, first achieved in 1985 by Merlin and collab- 
orators 0. In turn, the rich properties of these structures 
encouraged the study of systems based on alternative se- 
quences (e.g., Thue-Morse and hierarchical), which, in 
spite of not being quasiperiodic, still exhibit determinis- 
tic (or "controlled") disorder, i.e., they are nor random 
neither periodic. One of the most interesting features 
that many of these problems (either experimental or the- 
oretical) display is a fractal spectrum of excitations. For 
instance, experiments on Thue-Morse |3) and Fibonacci 
Q superlattices have uncovered scale invariant energy 
spectra. Numerical analysis of linear chains of harmonic 
oscillators with hierarchical nearest-neighbor couplings 
and equal masses exhibit spectra related to the triadic 
Cantor-set || (the Cantor-like structure is preserved even 
if masses are also distributed in a hierarchical way (^|). 
It has also been proved that the energy spectrum of a 
chain made of identical springs and of masses of two dif- 
ferent kinds arranged after the Thue-Morse sequence is a 
Cantor- like set J|. 

Even though controlled disorder typically leads to mul- 
tifractal spectra, in some cases only a few scales are suffi- 
cient for a satisfactory understanding of their thermody- 
namics. For instance, the phonon spectrum of the chain 
in [Q is essentially a one-scale Cantor-set; the electronic 
spectra that arise from Fibonacci tight-binding Hamilto- 
nians (either on-site or transfer) are governed by a couple 



of scale- factors |7| . 

Within this context, few-scale fractal spectra consti- 
tute simple prototypes for testing the thermodynamical 
implications of deterministic disorder. As a first step in 
this direction, in Refs. ^,|), one- and multi-scale frac- 
tal energy spectra were studied within Boltzmann statis- 
tics. It was shown that the scale invariance of the spec- 
trum has strong consequences on the thermodynamical 
quantities. In particular, the specific heat oscillates log- 
periodically as a function of the temperature. Moreover, 
general scaling arguments and a detailed analysis of the 
integrated density of states allowed for a quantitative 
prediction of the average value (which is related to the 
average density of states), period and amplitude of the 
oscillations. 

In this paper we extend the analysis of [|||] to TV- 
particle systems described by quantum statistics. Wc 
present results concerning the relevant thermodynamical 
quantities, i.e., chemical potential, average particle num- 
ber and specific heat. Our conclusion is that for phonons, 
and for bosons in general, the Boltzmann scenario sur- 
vives the inclusion of quantum symmetries. That is, 
the crudest (power-law) approximation to the density of 
states leads to a correct prediction of the average behav- 
ior of the specific heat. Log-periodic oscillations decorat- 
ing the average value are well reproduced when the sim- 
plest non-trivial corrections to the power-law density of 
states are considered. In the electronic case, the presence 
of a Fermi surface makes the thermodynamical quantities 
very sensitive to the particle number, thus preventing 
the use of smooth approximations. However, for special 
cases, interesting regularities can still be observed. 

The paper is organized as follows. In Section II we 
discuss some general aspects of the problem. We make 
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as well a short revision of previous results which will be 
useful in subsequent sections. In particular, we introduce 
the spectra that will be considered and describe their 
significant characteristics. In Section III we treat the 
phonon case, whose essential features can be explained 
analytically by means of scaling arguments similarly to 
the Boltzmann case. Section IV is devoted to Bose (gen- 
eral case) and Fermi statistics. Section V contains the 
concluding remarks. 



II. GENERAL CONSIDERATIONS 

We consider a spectrum constituted by levels lying on a 
{fi > T2) fractal set |l0| . Such spectrum can be constructed 
recursively starting from an arbitrary discrete set of lev- 
els {e' -'} belonging to the interval [0, A]. The hierar- 
chy {e( n+1 )} consists of two pieces which are obtained by 
compressing the set {e^} by factors r\ and r 2 , respec- 
tively. The lowest end of the first piece is positioned at 
e — 0, the upper end of the second piece is positioned 
at A (see Fig. 1). The fractal that arises in the limit 
n — > 00 has a (Hausdorff) dimension df given implicitly 

by r c [ s + r 2 f = 1- In the particular one-scale case, i.e., 
n = T2 = r, one has the explicit result df = — ln2/lnr. 
An alternative way of viewing the construction of the 
fractal is that in terms of gaps and "bands" . At hierar- 
chy n there are 2" — 1 gaps (in gray in Fig. 1), which define 
2™ "bands" (in white). When n is increased by one, the 
preexistent gaps remain untouched, but each band splits 
into two new ones and a new gap appears. 

We expect the thermodynamical properties associated 
to this spectrum to depend only on its hierarchical or- 
ganization and not on the specific pattern for n = ||. 
For this reason, and for the sake of notation simplicity, 
we will take {e (0) } = {0, A} (so that there are 2" +1 lev- 
els at generation n). In addition, we make A = 1 and 
fcs = 1 (ks stands for Boltzmann's constant) throughout 
the paper. 



A very useful tool in describing Boltzmann systems 
with fractal spectra |||| is the integrated density of 
states J\f(e) defined as 

JV(e) = f e de'D(e'), (1) 
Jo 

where D(e) is the one-particle state density. In spite of 
the very fragmented structure of the density of states 
(a set of delta functions with fractal support), a smooth 
approximation is usually good enough for computing in- 
tegral quantities, e.g., the specific heat. In the present 
case, the scaling properties of the spectrum imply that in 
a log-log plot Af(e) looks as a self-similar staircase with 
a fractional average slope. If first-order log-periodic cor- 
rections are included, we obtain a smooth approximation 
of the form Q 

N{e)/2 n+1 «e d [a + 6cos(wlne - (j))] . (2) 

The so-called spectral dimension d [jllj and the frequency 
uj arise from simple scaling arguments, yielding 

rf=-ln2/lnri ; w = -27r/mri • (3) 

The parameters a ~ 1, b, and <f> are functions of both 
T\ and r 2 and result from a detailed analysis of the ex- 
act spectral staircase. The smooth cummulative density 
allows for analytical manipulations while still keeping 
the principal ingredients of the self-similar spectrum. In 
fact, in the Boltzmann case, the use of (||) together with 
a perturbative approach leads to an accurate expression 
for the specific heat per particle as a function of the tem- 
perature T: 

c(T) K,d + a' cos(w In T - <f>) , (4) 

where a' <C d is a function of r\ and r 2 (9) . Eq. (|^) relates 
explicitly the log-periodic nature of the specific heat to 
that of the spectral staircase. 

III. PHONON STATISTICS 



n=0 



n=1 



n=2 



n=3 



1 







FIG. 1. The two-scale energy spectrum. Shown are the 
first hierarchies in the construction of a (ri,r-2) fractal spec- 
trum. The initial pattern (n = 0) was chosen as a set of two 
levels, {ei°',e2°'} = {A A}. At generation n, 2 n_1 new gaps 
(in gray) appear. The (positive) scale factors are restricted 
by the condition n + T2 < 1. 



A standard way to study the acoustic properties of a 
lattice is to consider a nearest-neighbor harmonic chain. 
This model is represented by the equation of motion 



m ; 



d Uj 
~aV~ 



(5) 



where Uj is the displacement of the jth atom (of mass rrij) 
from its equilibrium position, kj j±\ are the strengths of 
the couplings between neighboring atoms. Deterministic 
disorder is introduced by requiring the masses and/or the 
strengths to follow substitutional (e.g., Thue-Morse, Fi- 
bonacci, Fibonacci-class 12 1) or hierarchical sequences. 
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Assuming that the time dependence in (g) goes like 
Uj oc exp(— icot), the stationary equation of motion is ob- 
tained. Upon diagonalization one calculates the normal 
modes and the eigenfrequcncies u>. In the case of Thue- 
Morse ||, Fibonacci Jl3|— p~6^| and Fibonacci-class |l7]] se- 
quences the spectrum of energies e — hw has a Cantor- 
like structure only in the short wavelength regime. For 
long wavelengths the number of gaps and their sizes tend 
to zero and one recovers the standard results for the peri- 
odic lattices. On the other side, the hierarchical chains of 
Refs. ^||| do possess spectra with uniform scaling of the 
type (ri,r2). For these cases we can apply the formal- 
ism of J9| to obtain an explicit expression for the specific 
heat. The specific heat per atom c(T) is calculated as 



c(T) 



2™+' 

— y 

n+l Z-^i 



ej /2T 



sinh(ej /2T) 



(G) 



where the summation runs over all energy levels of a 
(n.) ^2) fractal of n-th generation. In Fig. 2 we illustrate 
the dependence of the specific heat (||) with temperature 
for different values of (ri,^). Logarithmic scales have 
been used to clearly display that, in the low temperature 
regime (T <C 1), a) c(T) increases in average as power 
law; b) there are log-periodic corrections to the average 
power law. As T is lowered this behavior persists down 



to a minimum temperature T n 



corresponding to 



the smallest scale of the fractal truncated at hierarchy n 
(when n — > 00, T m i n — > 0). Of course, for high temper- 
atures (T ^> 1), the classical result c(T) = 1 for a one 
dimensional lattice is recovered. Features a) and b) can 
be shown to be a direct consequence of the log-periodicity 
of the spectral staircase, which is captured by the smooth 
approximation (^). In order to show that this is the case 
let us first rewrite (0) as an integral: 



z{T) = 



1 



2 n+l 



e/2T 



sinh(e/2T) 



djV(e) 



(7) 



Now we make two approximations in Eq. (0). First we 
replace the exact j\f(e) by its smooth version (|J). Sec- 
ond: as our interest is the low-T regime, we extend the 
upper integration limit to infinity (this is equivalent to 
considering an unbounded fractal spectrum) . After some 
simple algebra we arrive at 



c(T) w T d d [a" + b" cos(w InT - , 
c" sin(u> In T — <f>)] . 



(8) 



Here a" = aa', b" = bib' - c'oj/d), and c" = -b(c' 
b'uj/d); the parameters a', 6', c! are given by 



dx ■ 



sinrr(x/2) 



1 

cos(wlna;) 
sin(w lnx) 



(9) 



Eq. (||) is the phonon analogue of the Boltzmann result 
(0) . It displays explicitly the average behavior of the spe- 
cific heat and the log-periodic corrections. In Fig. 2 we 
have plotted the analytical curves (||) for the same set 
of scaling factors used for the numerically exact calcula- 
tions. 
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FIG. 2. Phonons in a two-scale fractal spectrum. (a) 
Log-log plot of the specific heat per atom versus temperature. 
The curves are parametrized by the scale factors (r\,r2) in- 
dicated in the figure, (b) c(T)/T d vs InT. In all cases n = 9. 
Full lines correspond to numerical results (^) and dotted ones 
to our analytical prediction (|^), which is valid only for low 
temperatures (this is the reason why the analytical curves do 
not saturate at high temperatures). 

For low temperatures, the agreement between the ex- 
act and analytical results is very good [notice that the 
disagreement in the mean value for the (0.30,0.15) case 
actually corresponds to a 3% error]. The scale factor r\ 
is responsible for the period of the oscillations through 
uj = — 2tt/ lnri (as shown in Fig. 2b the period decreases 
as r\ increases). Moreover, as T\ controls the gap-scaling 
with respect to T — 0, smaller r\ values imply compara- 
tively bigger gaps, and, in turn, bigger amplitudes. For 
fixed T, decreasing r\ causes the spectrum to shrink to- 
wards e = 0, as a consequence the specific heat grows. 
These are the features observed by Petri and Ruocco 
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in their paper on one dimensional chains with hierarchi- 
cal couplings . We can now understand quantitatively 
their results as being a direct consequence of log-periodic 
corrections to the pure power-law scaling of the density 
of states (|[) . 



of high temperatures and small N, one obtains from ([Hj) 
z — > N* = N/2 n+1 . In Fig. 3b we have chosen a small 
value of r\ to show that [i indeed oscillates as the tem- 
perature gets through the scales of the spectrum. In this 
sense, Fig. 3b can be thought as an amplified version of 
Fig. 3a. 



IV. BOSE AND FERMI STATISTICS 

The statistical problem of particles with non identi- 
cally null chemical potential stands at the next level of 
complexity, both from analytical and numerical points of 
views. As in this case analytical derivations seem not fea- 
sible (except in limiting regimes), we resort to a numer- 
ical procedure. For a fixed average number of particles 
N we extracted the chemical potential /jasa function of 
the temperature by solving the equation 



N 



1 



2 „+i 

^2 e (ej-lA/T ± 1 



(10) 



where the sum runs over the levels of the fractal trun- 
cated at hierarchical depth n. The sign +/— corresponds 
to the Fermi/Bose case. (We are not taking into account 
spin degeneracies, which will not modify our results.) 
The numerical solution of Eq. ( |Io| ) is not a difficult task, 
however some care has to be taken with the choice of an 
appropriate initial value for /i. Once the chemical poten- 
tial has been obtained, the specific heat can be calculated 
as the T-derivative of the average total energy E(T), 



E(T) = J2 



(11) 
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A. Bosons 

Even though bosonic excitations with non-null chemi- 
cal potential do not seem to be of main relevance for the 
thermodynamics of superlattices, we decided to include a 
discussion on this topic for reasons of completeness, and 
because bosons in a fractal spectrum display a paradig- 
matic behavior. 

The dependence of the fugacity z = exp(/i/T) with 
temperature is illustrated in Fig. 3a. Several curves, 
parametrized by the average particle number (indicated 
in the figure) have been drawn. First of all, Fig. 3a dis- 
plays the following limiting features. As T decreases, and 
the system evolves towards condensation, the chemical 
potential increases from negative values up to a limiting 
value close to zero (we recall that our lowest energy is 
e = 0). The number of particles in the condensate is 
given by N$ = z/(l — z). Thus, for very low T, together 
with large N, one has z « 1 — 1/N . In the opposite limit 



FIG. 3. Bosons in a two-scale fractal spectrum. Fugacity z 
vs. InT, for different values of particle number TV (indicated 
in the figure). The chosen scales are (r^rs) = (1/3,0.2) in 
(a) and (0.02,0.2) in (b). In both cases n = 7. 

The dependence of the specific heat per particle on 
the temperature is illustrated in Fig. 4. Shown are the 
exact result and the numerical calculation which uses 
the smooth expression (||) for the spectral staircase. For 
low temperatures the agreement is excellent. It is clear 
that the mean value of the specific heat is indeed asso- 
ciated to the average level density and that it suffices 
to take into account only the first non-trivial correction 
to the power-law scaling to explain the oscillations. [It 
is worth pointing out that this scenario breaks down in 
the fermionic case (see later)]. We have tested that the 
curves in Fig. 4 do not change in the thermodynami- 
cal limit of increasing both N and the hierarchy n but 
keeping N* = N/2 n+1 fixed, except that the oscillatory 
behavior extends to lower temperatures of the order of 
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r" (the total number of oscillations of the specific heat 
is equal to the depth n). The relative particle number 
N* plays the role of a density, given that the number 
of levels below a certain energy grows with the volume 
of the system. When N* is sufficiently small and the 
temperature is high enough, the curves tend to those for 
the Boltzmann statistics This is analogous to the 

usual statement that the parameter X 3 p (where A and p 
are the thermal length and the density, respectively) tells 
how good is the classical approximation. 
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FIG. 4. Bosons in a two-scale fractal spectrum. Log-log 
plot of the exact specific heat per particle c(T) vs. T, for dif- 
ferent values of average particle number N (full lines) . Dotted 
lines represent a calculation which uses the cummulative den- 
sity of states (^). The horizontal line is c(T) = d, which is 
the classical average value for low T. At high temperature, 
the specific heat decays according to the classical 1/T 2 law. 
Scale factors are (ri,ra) = (1/5, 1/3). The fractal has depth 
n = 9. Base- 5 logarithms were used to show that there is one 
oscillation per "decade". 

Although an analytical description of the problem 
is lacking, our numerical calculations show that log- 
periodicity is robust enough to resist the inclusion of 
bosonic symmetries together with the restriction of par- 
ticle conservation. 



where ipj denotes the wave function at site j and {tj} 
are the hopping matrix elements. If these hopping ele- 
ments are arranged according to the Fibonacci sequence, 
the spectrum of energies {e} is essentially a fractal of the 
type (7*1, r%) (see e.g. 0), independently of the boundary 
conditions. However, the simple scaling is lost if the more 
general Fibonacci-class |l2| sequences are considered. 

The numerical scheme used for computing the ther- 
modynamical quantities in the bosonic case can be easily 
adapted to the electronic problem. Figs. 5 and 6 exhibit 
some relations among chemical potential p, average num- 
ber of particles N and temperature for the r± = r% = 1/3 
case. Fingerprints of the fractal spectrum are clearly seen 
in the dependence of p with temperature (Fig. 5). In 
the limit of zero temperature, for a given integer particle 
number N, p takes the value corresponding to the mid- 
dle of the gap (eAT,£Ar + i), as happens, e.g., in intrinsic 
semiconductors. In order to keep N fixed as tempera- 
ture grows, the Fermi surface moves in the direction of 
the lowest density of levels, integrated over an energy 
interval T. The relative concentration of levels in the 
neighborhood of the gap is a fluctuating function of the 
scale UbT. This gives rise to an oscillatory process that 
persists until T overcomes the largest gap (log 3 T « —1). 
From this point on, if N* > 0.5 (N* < 0.5), p tends in a 
monotonic way to +00 (—00). 



0.6 



0.4 



(J. 0.2 







-0.2 
























log 3 T 



B. Fermions 

The results presented in this section are valid for 
fermions in general, however, we assume that the frac- 
tal (ri, T2) corresponds to the electronic spectrum of a 
certain superlattice. From the theoretical point of view, 
such spectra appear in the simplest model for studying 
the electronic properties of a lattice, i.e., a stationary 
tight-binding equation, for instance, in its transfer ver- 
sion: 

tj+ii'j+i + tji'j-i = eipj , (12) 



FIG. 5. Fermions in a fractal spectrum. Chemical poten- 
tial /1 as a function of log 3 T for all integer particle numbers 
1 < N < 2 n+1 = 128 (increasing from bottom to top). We 
have only plotted the lower part of the full figure because 
it is symmetric with respect to p — 0.5. Parameters are 
ri = T2 = 1/3, n = 6. 

For smaller values of r± , oscillations in p are magnified 
(not shown), analogously to what happens with bosons. 

Fig. 6 displays the chemical potential as a function of 
the reduced number of particles N* for three values of 
fixed temperature. This figure gives some insight of the 
difficulties involved in solving Eq. (ITol). For low (fixed) 
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temperatures, fJ,(N) tends to the spectral staircase and 
the numerical problem gets relatively difficult. As tem- 
perature grows the staircase gets smoother and numerics 
simplify. We have taken N to be continuous, but observe 
that for integer N and low temperatures (compared to 
the gap size) the curves pass through the middle of the 
gaps, as it is clear in Fig. 5. 
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FIG. 6. Fermions in a fractal spectrum. Log-log plot of 
the chemical potential fi versus the reduced number of par- 
ticles N* for fixed temperature. The curves correspond to 
T — 3~ 9 , 3~ 6 , and 3" 3 (smoother curves correspond to higher 
temperature). Gray lines indicate the energy levels. We used 
n = ra = 1/3 and n = 4. 

The specific heat in the fermionic problem depends 
strongly on the position of the Fermi surface. Thus, we 
begin by analyzing some simple cases where the number 
of particles takes special values. For this purpose it is 
useful to recall the picture of bands and gaps of Fig. 1. 
Each gap determines a special number of particles N*, 
namely, the relative number of particles that would lie be- 
low that gap at T = 0. For instance, the first-generation 
gap determines N* = 1/2. The n — 2 gaps correspond 
with TV* = 1/4,3/4. Conversely, the denominator and 
the numerator in N* indicate respectively the generation 
n of a gap and the number of filled bands of nth genera- 
tion, e.g., N* = 1/32 = 1/2 5 corresponds to a gap which 
appeared in the fifth generation and to one filled band of 
fifth generation. 

A first example that illustrates the structure of the spe- 
cific heat is depicted in Fig. 7. There we consider a spec- 
trum with two different scales r\ = 1/3, r2 = 1/9, and the 
selected set of numbers N* = 1/32,1/16,15/16,31/32. 
For the first two cases, N* = 1/32,1/16, the size of 
the gap is of the order of the Fermi temperature Tp. 
Once the temperature overcomes the width of that gap, 
the electrons can fully access the next band above, so 
that the mean level occupation drops approximately to 
1/2. From then on the electrons behave essentially as 



Boltzmann particles: the specific heat oscillates around 
a constant average value, as expressed by Eq. (Q). The 
average value and the period are governed by the scal- 
ing factor n = 1/3 through Eq. (0). The number of 
oscillations is equal to the generation index n: four if 
N* = 1/16 and five if N* = 1/32. The case of numbers 
N* = 31/32, 15/16 is equivalent to having respectively 
N* = 1/32, 1/16 holes and it is then complementary to 
the previous one. Remarkably, in the case of holes, the 
dominant scale is r-i. The reason for this is that ri rules 
the scaling with respect to the upper edge of the spec- 
trum, and, in fact, with respect to the upper (lower) 
edge of every band (gap); in other words, r-i is the rele- 
vant factor for negative temperatures. Eqs. fj|) and (^) 
are still valid provided that r\ and are interchanged. 
According to these arguments, the choice r-i = r\ implies 
that the period of the oscillations for holes is twice that 
for electrons, which is confirmed by our calculations (see 
Fig. 7). Horizontal lines correspond to the Boltzmann 
average result, properly normalized in the case of holes. 
Increasing the hierarchy n and keeping N* fixed does not 
change Fig. 7 because the gap on top of the filled band 
determines the smallest scale of the fractal that can be 
resolved. 
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FIG. 7. Fermions in a fractal spectrum. Log-log plot of 
the specific heat versus temperature. The reduced numbers 
of particles are iV* = 1/32, 1/16, 15/16, 31/32, and the scale 
factors ri = 1/3, r2 = 1/9. The horizontal lines represent the 
classical average value for the appropriate number of "hot" 
electrons (N* = 1/32, 1/16) or holes (N* = 15/16,31/32). 

Other special particle numbers allow for observing a 
mixed electron-hole behavior. Let us consider, for in- 
stance, the cases N* — 7/64,57/64, corresponding, re- 
spectively, to 7 and 57 filled bands of sixth generation 
(Fig. 8). For N* — 7/64, the Fermi temperature is several 
scales bigger than the width of the seventh gap, so that 
the low-T part of the specific heat is associated to hole 
excitations jumping from the (empty) eighth band down 
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over the sixth-, fifth-, and fourth-generation gaps. These 
excitations are of Boltzmann nature and govern the spe- 
cific heat until the temperature becomes of the order of 
Tp. At this point electrons are able to jump up over 
the third-generation gap between the eighth and ninth 
bands. The high-T (T > 3~ 4 ) oscillations are associated 
to electronic excitations through the third-, second-, and 
first-generation gaps. Similar considerations apply to the 
complementary case N* — 57/64. The horizontal lines in 
Fig. 8 correspond to our predictions for the average spe- 
cific heat, which take into account the effective number of 
electrons or holes that contribute to the specific heat in 
each regime. (The fact that the case N* — 7/64 presents 
one less oscillation than its complementary is due to an 
overlap of temperature scales in the transition from hole 
to electron behavior.) 
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FIG. 9. Fermions in a fractal spectrum. Log-log plot of the 
specific heat versus temperature for the numbers N = 25, 77, 
and the scale factors n = 1/3, ri = 1/9 in = 7). 
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FIG. 8. Analogous to Fig. 7 but for the particle numbers 
N* = 7/64, 57/64 (see text). 



The general case of an arbitrary number of particles 
shows a non-trivial mixture of the simple behaviors de- 
scribed before. However, the gross features of the specific 
heat for arbitrary iV can be qualitatively (and sometimes 
quantitatively) understood as being a reflection of se- 
quences of electronic and hole excitations that alternate 
themselves in producing the oscillatory patterns shown 
in Fig. 9. 



Formulations based on smooth approximations to the 
density of states, analogous to those made for bosons, are 
bound to fail in the fermionic case. The intrinsic discon- 
tinuity of the Fermi problem, which is critically enhanced 
by the dual scaling respect to the lower and upper edges 
of the gaps, is in essence non compatible with smooth 
approximations (with the exception of too special cases) . 



V. CONCLUDING REMARKS 

We have analyzed the quantum statistics of model 
systems exhibiting two-scale fractal spectra, with spe- 
cial emphasis on the structure of the specific heat. Our 
findings extend previous results on classical statistics to 
show that the thcrmodynamical manifestations of spec- 
tral fractality are robust with respect to the inclusion 
of quantum symmetries. The general scenario for Boltz- 
mann and bosonic particles can be summarized as fol- 
lows. In spite of the very fragmented structure of the 
real density of states, a formulation which starts from a 
smooth approximation but takes into account the coars- 
est log-periodic fluctuations is sufficient for a good de- 
scription of the specific heat and other averaged quanti- 
ties. In the case of fermions, even if a treatment which 
is uniform in the particle number is not possible, many 
features of the problem can be understood with the help 
of the simple results for the Boltzmann case. 

Our analysis was limited to two-scale fractals. Previ- 
ous results H indicate that provided the scaling towards 
the inferior limit of the spectrum is uniform, the inclu- 
sion of additional scales will not change our conclusions 
in what concerns the bosonic case. Fermions, however, 
can feel the scaling with respect to any point of the spec- 
trum. So, the complexity of the fermionic problem would 
be proportional to the number of relevant scales, in spite 
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of uniform zero-energy scaling. 

Let us finish by mentioning that the log-periodicities 
described in this paper may be observable in real phys- 
ical systems, e.g. in the Fibonacci superlattices p|,|l8|]. 
In fact, although our discussion was restricted to perfect 
deterministic systems, both experiments and numerical 
simulations by Todd et al. indicate that the hierar- 
chical organization of the electronic bands may be pre- 
served even if substantial amounts of (random) disorder 
are added to the system. From a more general point of 
view, Saleur and Sornette (T^| have demonstrated that 
the connection between discrete scale invariance and log- 
periodic oscillations is robust with respect to the presence 
of disorder. 
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